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Abstract. 

This paper is concerned with the stability of periodic wave trains in a gener- 
alized Kuramoto-Sivashinski (gKS) equation. This equation is useful to de- 
scribe the weak instability of low frequency perturbations for thin film flows 
down an inclined ramp. We provide a set of equations, namely Whitham's 
modulation equations, that determines the behaviour of low frequency per- 
turbations of periodic wave trains. As a byproduct, we relate the spectral 
stability in the small wavenumber regime to properties of the modulation 
equations. This stability is always critical since is a 0-Floquet number 
eigenvalue associated to translational invariance. 
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1 Introduction 

We study perturbations of periodic travelling wave solutions to the following 
generalized Kuramoto-Sivashinski (gKS) equations 



where u{x,t) G M is a heigth at place x G M and time t > 0. Actually 
this equation describes approximately long and small surface waves of two 
dimensional incompressible viscous fluid down an inclined plane. Here 6R 
{6 > and R > 0) measures the deviation of the Reynolds number from the 
critical Reynolds number above which long wave perturbations are spectrally 
unstable. Since R > 0, steady flows are unstable and periodic travelling 
waves, so called roll-waves, appear. 

This transition to instability is proved to occur also in shallow water 
equations (SW): 



(here h{t, x) G M is an heigth, u{t, x) G M a velocity, F > a Froude number 
and -Re > a Reynolds number). If F > 2, the stationnary solutions are 
unstable and roll-waves also appear. These roll-waves are observed both 
experimentally and numerically. 

Both in ([1]) and ([2j), roll- waves are proved to exist through a Hopf bi- 
furcation argument and a whole family of periodic solutions exists, which 
ends up with a solitary wave. This family is parametrized by two quanti- 
ties, period and either spatial mean heigth for ([T]) or relative discharge rate 
for ([2]). Under some conditions, a similar situation is also proved to ex- 
ist in the Navier-Stokes equations (NS) for an incompressible fluid down a 
ramp: indeed, taking into account capillarity and under a spectral assump- 
tion, namely a Hopf bifurcation scenario, it is proved that small amplitude 
periodic travelling wave solutions to (NS) exist when R^. > 5 cotan(6')/6, 6 
being the angle of the slope [TT]. 

Actually the (gKS) equation is a generic equation which describes this 
transition to instability in the regime of small amplitude long waves. One can 
derive formally (gKS) either from (SW) or (NS): see e.g. [16] for a derivation 
from (NS) and [T7] for a derivation from (SW). These formal derivations 
are strongly connected to the issue of the rigorous derivation of the shallow 
water equations from the Navier Stokes equations in the case of viscous fluids 
flowing down a ramp. Assuming capillarity, the (SW) equations are derived 
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rigorously from (NS) whenever steady solutions are stable The (gKS) 
equations is the reference equation to describe small amplitude solutions of 
(NS) or (SW) when the flow is weakly unstable. In the stable regime, both the 
shallow water equations and (NS) can be approximately reduced to a single 
viscous Burgers equation (or Benney equation): see [15] for the reduction 
from (NS), as for (SW) the limit is a standard relaxation limit in balance 
laws. Somehow this reduction already provides us with a validation of (SW) 
in the stable regime. Yet solutions to the Burgers equation blow up when Re 
exceeds a critical value and (gKS) is then needed. 

Stability of roll-waves and description of the dynamic around such wave 
trains are interesting intricated questions. Some information may already be 
deduced from the fact that roll- waves emerge from unstable steady solutions. 
At the onset of roll-waves, that is in the small amplitude, the wave train is 
close to a stationnary solution and so is expected to be spectrally unstable. 
For the same reason, the solitary wave is also expected to be spectrally 
unstable, so that only a range of periodic roll-waves could be stable. This 
would explain the coalescence dynamic that is usually observed (see e.g. [T]). 
To go beyond this general heuristic discussion, we summarize here some of 
the results obtained for the (SW) equations in [2], [3l H] . It is always a non 
trivial task to carry out the spectral analysis of equations linearized about 
periodic travelling waves. Yet, using perturbation methods or reducing to 
numerical computations in a finite box through energy estimates or stable- 
unstable tracking, solitary waves and small amplitude wave trains are there 
proved to be spectrally unstable and a range of periodic roll-waves is indeed 
proved to be stable. The low frequency instability is a strong instabihty 
for small amplitude (the spectrum crosses transversally the imaginary axis), 
whereas, for solitary waves, essential spectrum is unstable (corresponding to 
steady states spectrum) and point spectrum is weakly unstable (the spectrum 
is tangent to the imaginary axis but on the unstable side). Note that the 
spectral stability of some periodic wave trains may be used to prove also 
their nonlinear asymptotic stability [lOJ. 

Complementary to direct numerical computations, another approach to 
study the stability of periodic wave trains is to carry out a long wavelength 
analysis in the neighbourhood of roll-waves since this is the particular regime 
governing the return to equilibrium. Slow modulation of roll-waves are then 
proved to be (still) described locally by roll-waves whose parameters evolve 
on a slow scale according to a set of first order partial differential equations, 
so called Whitham's averaged equations. This approach involves at least two 
steps, the derivation of averaged equations and the proof of a connection 
between spectral proprerties of both the modulation system and the original 
equations. Obviously the spectral step yields then some necessary conditions. 
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expressed through averaged quantities, for the stabihty of periodic travelhng 
waves. Such an approach was proposed for conservation laws in [H] where 
it is proved that the leading term, in the expansion about the origin, of the 
corresponding Evans function is given by the dispersion relation associated 
to linearized Whitham's equations. For (SW), an Evans function approach 
to stability was also introduced in [12] and extensively used in [21 [21 II], and 
his relation with the dispersion of the Whitham's system also established in 
[13] . Yet, in [13] a Bloch transform approach is shown to give a more nat- 
ural and powerful way of achieving the spectral step. Actually, in [13], the 
authors obtained, as in [13] for conservation laws, that the tangency to the 
imaginary axis of low frequency spectral curves is descibed by the hyperbol- 
icity of a first order modulation system but also that the curvatures of this 
curves are related to the parabolicity of a second order modulation system 
also derived in [13]. Moreover, the Bloch transform approach gives almost 
readily a relation not only between eigenvalues but also eigenvectors. This 
is precisely the keystep in order to perform a third step in the modulation 
analysis : to validate at a nonlinear level the modulation equations as provid- 
ing a good approximation of low-frequency perturbations of wave trains as 
is explained in [6j for reaction-diffusion systems (in cases where modulation 
yields a single scalar equation and not a system). In [13] the authors proved 
such a nonlinear validation of the first order modulation system when its 
hyperbolicity is met. An expansion of eigenvectors would certainly also be 
useful for a justification of second order modulation systems and construc- 
tions of generalized shocks solutions corresponding to modulated roll-waves 
with shocks in parameters as performed for scalar modulation equations in 
reaction-diffusion systems [6]. 

In this paper, we carry out the derivation and spectral steps of a similar 
analysis for (gKS). For fixed 6 > 0, the set of periodic traveling waves is 
two dimensional (up to translations) and it is proved to exist through a 
Hopf bifurcation argument. This set is parametrized by the wavenumber 
and spatial mean. In the low frequency regime, we compute a set of two 
equations that governs slow modulation of the local wavenumber and local 
spatial mean. Similarly to roll- waves in (SW), we obtain "inviscid" (first 
order) and "viscous" (second-order) models and relate these systems to the 
spectral stability of periodic wave trains just as in \13\ . 

However, recall form the derivation of (gKS) from (NS) and (SW) that 
the parameter 6 should be small so that the domain of validity of Whitham's 
equations, whose size is measured as the size of the allowed perturbation 
characteristic wavenumber in the modulation, shrinks to 0. One has to carry 
out a new modulation analysis with characteristic wavenumbers e comparable 
to 6. We do this setting 6 = 6e and we show in this case that the Whitham 
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set of equations is composed of three first order partial differential equations 
with a balance term which takes into account the dissipation of (gKS). The 
number of equations, three, is easily understood since, in the limit 5 — )■ 0, 
(gKS) is a perturbation of the Korteweg-de Vries (KdV) equation which 
possesses a three dimensional set of periodic wave trains. We will connect 
this set of modulation equations with the spectral stability of periodic wave 
trains in the regime of wavenumber of order 6. We also relate the modulation 
system with the one for (KdV) (in the limt 5—7-0) and the one obtained for 
a fixed S (in the limit 5 — )■ cxd). 

The paper is organized as follows. In section [2], we recall the structure 
of periodic wavetrains solutions to (gKS) both in the case 6 > fixed and 
in the limit 5 — )■ 0. In section [31 we consider the case 6 > fixed, we 
first derive "inviscid" and "viscous" modulation equations, then carry out 
a direct spectral analysis of (gKS) and relate it with a spectral analysis of 
the modulation systems. In section HI we consider the case 5 — )■ 0, compute 
the set of Whitham's equations, carry out a direct spectral analysis in the 
regime of frequencies of order 5 — )■ and again show the relation with the 
Whitham's systems. 

2 Periodic travelling waves in (gKS) 

Up to changes 

t = R^/h, X = R^/^x, u = Ru,6 = R^''^5 

one may assum^ R = 1 and consider (gKS) equation in the form 

dtu + Qud^u + dlu + 5 {dlu + d^u) = . (3) 

For the moment we search for periodic travelling waves, with profile U a 
periodic function, without fixing the period of the profile, and thus for u in 
the form u{x, t) = U{x — ct). This yields 

- cu' + mu' + U" + 5 [U" + U"") = . (4) 

Integrating once, we search for a periodic profile U, a speed c and a constant 
q such that 

-cU + 3U^ + U" + 5{U' + U"') = q. (5) 

For arbitrary 6 > 0, the existence of periodic solutions to @ may be 
obtained through a Hopf bifurcation argument. Equation ([5]) possesses two 

'''One could also eliminate 6 by rescaling u again but we keep it to fit with equations 
as ■written in [I]. 
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stationnay solutions t/_(c, q) < ?7+(c, q) such that 3U^—cU± = q. Linearizing 
at f/ = f/± yields 

{6U±-c)U + U" + 6{U' + U"') =0. (6) 

The assocatied characteristic equation is 5(A^ + A) + + (6f/± — c) = 0. 
It is easily seen that the stationnary solution undergoes a Hopf bifurcation 
whenever 6U± — c = 1. This can only happens for U = f/+ and the equation 
settles a CHopfiq)- Denote V = U — U+ and 6[/_|_ — c = 1 + /i with -C 1. 
Equation ([5]) then reads 

l^(3\/ + 1 + /i) + l^" + 6{V' + V") = . 

For fi sufficiently small on one side of 0, there is a family V{fi) of amplitude 
0{>/\Jj]) and frequency k{n), with k{0) = 1 and fc'(O) 7^ 0. 

As a result, one obtains a two dimensional manifold of periodic traveling 
waves (identified when coinciding up to translation) parametrized by the 
wavenumber k and q. Note that in the small amplitude regime one may 
also parametrize by k and the spatial mean of the solution. Indeed, there 
stands U =< U > +0{^/\Ji\), thus since q = 3 < > —c<U>, one 
finds 3 < U -c < U >= q + C(|yu|). At last recall that, for /i small 
enough, 6 < f/ > — c 7^ so that one can switch from q to M =< U > in 
the neighbourhood of the Hopf bifurcation. As a result one may prove the 
following proposition. 

Proposition Let M G M &e fixed and k < 1 be such that 1 — k is small. 
Then there exist a unique c{k, M) and a unique q{k, M) such that there exist 
a 1-periodic U{- ;k, M) solution to 

k (3f/2 - cU) + k^U" + 6 [PW + k^U'") = q, <U>=M . (7) 

Moreover this solution is unique up to translation. 

Maybe the last thing we should say about the proof of this proposition is 
why the bifurcation occurs for k < 1 (rather than A; > 1). A quick way to see 
this is to multiply the first equation in ([7]) by U' and integrate over a period. 
This yields < ([/')^ >= k^ < {U"Y > thus A; < 1. Note also that above and 
through the text we use < / > to denote the average of / over one period ; 
for 1-periodic functions this is just /g /. 

Let us now explain how to reduce the search of solutions to ([7]) to the 
case M = 0. Inserting the ansatz 

U = M + [/(°) (8) 
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into yields < U^^^ >= and again 
with 

c(°) = C-6M, = q + cM-3M^. (9) 

Note that in the {k, M) parametrization of profiles, none of the above quan- 
tities c*-°\ U^'^\ g*-^-' depends on M. From now on all quantities with a (0) in 
superscript will refer to the zero mean problem. 

For applications, we usually assume that 6 is small. Then equation 
is a singular perturbation of (KdV) equation whose periodic solutions are 
described with the help of elliptic functions. In the singular limit 5 — )■ 0, 
periodic solutions to (jl]) were proved in [7j to be close to periodic solutions 
to (KdV) and an expansion with respect to S was obtained. Moreover, in [T] , 
a formal spectral analysis was carried out, but only in the context of pertur- 
bations with zero mean. This does induce a lack of generality. Here, we will 
provide a spectral stability analysis for arbitrary (bounded) perturbations. 

We provide now expansions of roll-waves profiles in the 6-^0 limit. As 
5 — )■ 0, one may expand c(°), f/(°), q^^^ as 

ci^) = co + 6ci + Oi6^), = go + % + C'(5'), U^""^ =Uo + 6Ui + 0{6^). 
The 1-periodic solution to ([5]) for 5 = and < f/o >= is given by 

Co = 16Kipfne(2-p'-3§^), (10) 
go = 4(2ir(p)A:)^(-3|M^ + 2(2-p2)fg+p2-l), 

with p g]0, 1[, K{p) and E{p) elliptic integrals of first kind and second kind 
and dn( ■ ) Jacobi's delta function with modulus p. Up to this order, there 
is no selection of a particular wave train {p is arbitrary) and the manifold of 
periodic solutions (identified when coinciding up to translation) would be of 
dimension 2, whereas it was 1-dimensional for 6 > fixed when restricting 
to zero mean solutions. To recover the same dimension, one has to compute 
an expansion of the solution to the next order: one finds 

k^Lf^' + {6UoUi - doUiY - ciU'^ + klJ'^ + k^'U'^" = . (11) 

d ^ ^ 

The linear operator / = k"^—-^ + -r ((S^^o — cq) ■) is Fredholm of index 

dx'^ dx 

and 1 and Uq span the kernel of its adjoint so that one can readily deduce 
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that equation flTT]) has a solution provided that the following compatibility 
condition is satisfied, < (f/g)^ >= < (^o)^ >• This equation yields the 
selection criterion 

^ 16^(p) (12) 



where 



A;2 

6N{p) 



N{p) = 7(1 - p2 + p'^)E{pf - (10 - 15p2 + I3p4 - Ap^)E{p)K{p) 



7Dip) 

i„4 ,(„6\ TT-Z^N r^/^N 

(13) 

+ (3 - 6p2 + - 2p^)K{pf , 
L'(p) = 2(1 - ^2 + p^)E{p)/K{p) - (2 - 3j92 + p4) . 



In order to determine Ci, one has to consider higher order corrections to U^^'^: 
in fact, Ci is determined through a solvability condition on the equation 
for U2. This yields ci = (see [7] for more details). As a consequence, 
coming back to the more general case where M may not vanish, we have 
also obtained a two dimensional manifold of (asymptotic) periodic solutions 
(identified when coinciding up to translation) parametrized by spatial mean 
M and wavenumber k (or alternatively the parameter p g]0, 1[). Note that 
the limit /c — )■ (i.e. p — )■ 1) corresponds to a soliton and k 1 (i.e. p — )■ 0) 
corresponds to small amplitude solutions (or equivalently to the onset of the 
Hopf bifurcation branch) . 

In what follows, we establish Whitham's equations that describes modu- 
lation of these wave trains through slow evolution of their parameters. For 
fixed S > 0, this will give us a system of two partial differential equations. 
To analyze the limit 5 — > and obtain a consistant Whitham's modulation 
theory, one has to do something slightly different and choose a characteristic 
perturbation wavenumber e of order 0{6). The equation (gKS) is then a per- 
turbation of (KdV) and we obtain a set of three partial differential equations 
which are balance laws but with a source term. As — ?■ 00, this source 
term is stiff and the system relaxes to a system of two partial differential 
equations similar to the one written for fixed 6 > 0. 

Let us also introduce the notation Q = —kc for the rest of this work. 
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3 Stability of wave trains for fixed ^ > 

3.1 Modulation equations 

3.1.1 First order Whitham's equations 

Let us denote e the characteristic wavenumber of perturbations. To study low 
frequency perturbations on scales {x/e,t/e) and obtain an averaged modu- 
lated behaviour, we first rescal variables as {X,T) = {ex,et). Then equation 
([5]) turns into 

dTU + 6udxu + e^dj^u + 5{edlu + e^dj^u) = 0. (14) 

In what follows, we will assume 6 > fixed (but small). We search for 
solutions to ([Hj) in the (formal) form 

u{X,T) = f^s'uk (^^^^;X,t\ with y,(X,T) = f^^'^ji^^T) , 

fc=0 ^ ^ j=0 

where Ui{y; X,T) are 1 periodic functions in y. Inserting this ansatz into 
ffT^ and setting Qq = 9t0o and = dx<t>o, one finds at the 0{e~^) order 

Vt^dyUQ + GkoUodyUo + k^d^uo + 6 {kldyUo + k^dyUo) = . (15) 

Denoting Mq =< Uq >, equation f[T^ implies Qq = Q{ko,Mo) and (up to 
some translation) 

Uo{y;X,T) = U{ko{X,T), Mo{X,T)){y) . 

Compatibility condition OtOx'-Pq = dxdx^Po yields the first equation of a 
Whitham's system 

dTko-dxmko,Mo)) = . (16) 
Recalling Q = ~kc, this equation is also written as 

dTko + dxikociko,Mo)) = , (17) 

Now we identify 0{e'^) terms. Setting = dx^Pi and ki = dx^Pi, one finds 

Co Ml + druo + dxiSul) + 3 {k^d^yyUo + kodxkod'^uo) 

+ S {2kodxyUo + dxkodyUo + Ak^dj^yyyUo + Qk^dxkod^uo) (18) 
+ n.dyuo + ki {dy{?>ul) + ?>kld'yUo + 5{2kodluo + ^kldluo)) = 
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where Cq denotes the hnear differential operator defined by 
[>Co/](l/; X, T) = [Cko(x,T),Mo{x,T)f{ ■ ; X, T)] (y) , 

Cu,Mf = mU{k, M) - c{k, M)) f) + A:^^ + 5 (k'^ + k'^j , 

C being thus the operator associated to the hnearization about U {k, M) of 
the profile equation in ([7]) differentiated once. The kernel of £ = Ck,M is one 
dimensional and spanned by U'{k, M). Indeed, is a semi-simple eigenvalue 
of £ of multiplicity 2 and dMU{k, M) = 1 is a generalized eigenvector since 
C[l\ = 6kU'. The adjoint of C, is Fredholm of index with a one 
dimensional kernel! spanned by 1. As a result, equation f|T8l) is solvable if 

dr <uo> +dx <3ul> = 0. (19) 

Note that using the notations (Q corresponding to the decomposition ([8]), 
one has 

c = 6Mo + c(°)(^o) and < 3ul > = q + cMq = q^'^\ko) + . 

This yields the following Whitham's system 

drko + dx {ko c{ko, Mq)) = 0, drMo + dx {q^'\ko) + 3M^) = 0. (20) 

This is the classical Whitham's system in the case of (gKS). Later on, we 
will relate the hyperbolicity of this system with a stability index that is used 
to give necessary conditions for the spectral stability of periodic wave trains. 

3.1.2 Second order Whitham's equations 

We now proceed with the modulation expansion and find an evolution system 
for ki and Mi =< ui >. This system will be useful to complete the spectral 
stability analysis of periodic wave trains. We first rewrite equation (|T8|) in a 
more convenient way. Differentiating ([7j) with respect to k and M yields 

C[dkU] = -dk^ U' - ((3f/2y ^ 3^2f;/./ ^ ^^2kU" + APU"")) , 

CIOmU] = -dnn U', <dkU>=0, < OmU >= 1 . ^^^^ 

As a consequence, denoting by a subscript the evaluation in {ko, Mq) (so 
that for instance Uq = Uq), equation f[T5]) may be turned into 

Co[ui - kidkUo] + i^i - kidk^) dyUo + OtUq + dxiSul) + Rq (22) 

^The fact that 1 lies in the kernel is directly related to the fact that (gKS) is a conser- 
vation law. 
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with Rq defined as 

i?o = 3 {k^d^yyUo + koOxkod^uo) + 6 {2kodxyUo + dxkodyUo) ^^^^ 

+6 [Ak^dj^yyyUo + 6 ^ A: Q ^ Q ) • 

Note that Rq depends only on (^0,^0) and has zero mean. From f[T^ . we 
know that there is a unique /o satisfying £o[/o] = dxUo + (9x(3mq) + Rq and 
< /o, f^o 0- Equation (122|) now reads 

Co[ui - hdkUo + /o- < /o >] + (^^1 - kidk^ + 6A;o < /o >) = . (24) 
Let us denote by Mi the function of (X, T) satisfying 

- kidk^o - MidM^o + 6ko< fo> = . (25) 
Then equation is equivalent to 

m = MiaM?7o + kidkUo+ <fo> -fo + ^if/o (26) 

for some (independent of y). Note that Mi appears to be Mi =< ui >. 
We now write a set of evoultion equations for (/ci,Mi). We first obtain an 
equation that governs the time evolution of ki by using the compatibility 
condition drOxfi = dxdTfi and equation fl25l) . 

dTki-dx[ki— + Mi^^ = -dx<6kofo>. (27) 

In order to obtain an equation that governs the evolution of Mi, we have to 
consider 0{e) terms in f|T4l) . 

drui + dx {Quo ui) + 6dxxUo = dy{---) . (28) 

The exact form of the right hand side in fl28|) is not needed here to com- 
pute the set of equations for {ki,Mi). Indeed, equation fl28l) imposes the 
compatibility condition 

dr <ui> +dx < Quo ui > +6dxx < uo > = . (29) 

Inserting (l26i) into (l29i) yields 

drMi + 9x(d[<3f/2>][A;i,Mi]) 

= dx<Quoifo-<fo>)>-5dxxMo 

where d yields differential with respect to parameters {k,M). This almost 
completes the derivation of a viscous Whitham's system. Indeed, let us 
introduce k = ko + eki and M = Mo + eMi. Next, consider (|T7|) +£: (l27l) on the 
one hand and (fT9l) +£ (!30|) . Then, neglecting C(e^) terms, these equations are 
written as a system of viscous conservation laws with 0{e) viscous terms. We 
shall relate these viscous equations with the curvature of the spectral curves 
associated to the stability of periodic wave trains in the small wavenumber 
regime. 
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3.2 Spectral stability of periodic wave trains 
3.2.1 Bloch analysis of spectral stability 

In this section, we analyse the spectral stability of a periodic wave train 
U{- ; /c*, M^) Denoting by a star subscript the evaluation in {k^,M^), in an 
adapted co-moving frame equation linearized about f/^ reads 

Otu + kA ((6f/^ - c^)u) + kldlu + 5 {kldlu + k^dlu) = 0. (31) 

This is shortly written as dtu + C^^u where £^ = Ck^^M^,. We carry out a 
Bloch analysis of this problem and search solutions to (1ST]) in form u{x, t) = 
e^^~^'^^Uy{x), with ill, a 1-periodic function and u G zM. For writting conve- 
nience, we will drop the u subscript indicating u dependence. As announced 
we consider 

d W ..f d ^ ^^^^ 



+ 6 



k:\ ^ + 1^] u + kU — + u 
ax I \dx 



-Am . 



The value A = is an eigenvalue of £^ associated to the translation mode 
u = U' thus corresponding to z/ = 0. We compute the spectrum in the regime 
of large wavelength perturbations |z/| <^ 1 and expand solutions to equation 



([32]) as 

A = z/Ao + z/^Ai + C>(z/^), u = Mo + Z/Mi + Z/^M2 + C(z/^). 

We first identify O^u^) terms in equation fl32]) and find £*mo = so that 
Mo = koU^/ki, for some constant k^. Note that mo is mean free. 
Next, we identify 0(z^^) terms and obtain 

A«i + ^f/; + ko ((6f/, - cM + ?,kiu':' + 5 {2k,u': + Akiu:")) = o 

thus, using ( |2T]) and denoting Mq the constant satisfying 
A k 

+ kAc^ ko + kAic^ Mo = , (33) 

k^ 

we are left with 

£4mi - kodkU, - MoOmU,] = . (34) 
But equation flM]) is equivalent to 

Ml = kodkU, + ModMU, + ^Ul (35) 
12 



for some constant ki. Note that Mq =< ui >. In order to obtain a full set 
of equations for (A;o,Mo), one has to consider O^u"^) terms. We find that U2 
should satisfy 



A[m2] + AoMi + /c^ ( (6t4 - c^)ui + 'ik 



+ XiUq + okl— h k^uo + 6k 



dx dx'^ 







(36) 



dtjj y dijC 

Equation fl36|l has a solution if and only if 

Ao < Ml > +k^, < (6f/^ — c^)mi > = . 
Substituting ([35]) into this equation yields the equation 

hdk[< 3f/2 >], A;o + (Ao - he, + hdM[< 3f/' >].) Mq = . (37) 
Setting Ao = A;^Ao, the linear system fl33|37p has a solution if and only if 



afe[< 3f/2 >], Ao - + dM[< 3f/' >] 



. 



(3^ 



As a result, one obtains a necessary condition for the stability of periodic wave 
trains. Indeed, recall that A expands as A = vki.XQ + O^u'^) and z/ G iM so that 
for the periodic wave train to be stable, it is necessary that Aq G M. In next 
section, we will find out this former condition also expresses as hyperbolicity 
of the first order Whitham's system. 

Again we proceed with the computation of the expansion of eigenvalues. 
Let us first introduce and f)*^ defined as 

v'' = (Ao - c.) dkU, + 3 dk (f/2- < f/2 >), + 3 k,{U, + k.dkU,)" 

+5 (t/; + 6 klu': + 2 K{dku:)' + 4 klidku:)'") , 

v""' = (Ao - OmU, + 3 Dm {U^- < >), ■ 

Notice that v'' and v'^^ are mean free and therefore there are (unique) mean 
free and such that 

= k^v'' and C,[f''] = k,v'' . 

Then we denote Mi the constant (with respect to y) so that 

AiA;o + AoA;i + A;,dc,(A;i, M,) = -6k, (ko < > +Mo < f'' >) . (39) 
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Substituting fl5^ into equation fl5^ and using the compatibility condition 
flTr|) and relations ([5T]) yield then 

U2 = dU^ih, Ml) + ko{< f' > -h + Mo(< > + ^f/: (40) 

for some /c2 independent of ?/. Notice again Mi =< {^2 >. To complete the 
system on (A;i,Mi), we still have to consider (9(z/^) terms in (132|) 

C^us] + X0U2 + h ((6f/. - c^)u2 + Skl^ + 6 {2k f-^ + ^kf^^ ^ 

+ Aiui + kl {?,kf-^ + 5[u, + Qkl ^) ) (41) 

+ A2Uo + A;3uo + 45fc,^^ = 0. 

Setting Ai = Ai/fc^, the solvability condition for (HT]) expresses as 

(Ao - c*) < U2 > +Ai < Ml > + < QU^U2 > +Sk^ <u^> = . (42) 

Substituting ([35]) and (gO]) into (gS]) yields 

(Ao-cjMi + AiMo + 3d[<f/2>],(A;i,Mi) = -6k,Mo (43) 
+ < 6f/.(f - < /'^ >) > ko+ < 6f/.(/*^- < >) > Mo . 

It is easily seen by induction that we can obtain similary a complete set 
of equations for (A„,/c„,M„ =< itn+i >) for any n e N and thus obtain an 
expansion of both eigenvalues and eigenvectors up to any order in the small 
wavenumber limit z/ — > 0. We relate the two mode problems fl33f37l) and 
fl39f43l) with respectively "inviscid" and "viscous" Whitham's equations in 
the following section. 

3.2.2 Spectral validity of Whitham's equations 

In this section, we show the connection between the spectral analysis carried 
out previously and the spectral analysis of linearized Whitham's equations. 
Let us first check the validity of "inviscid" Whitham's equations, 

dTk + dx{kc{k,M)) = 0, dtM + dx {< SU"^ >) = 0. (44) 

Linearizing (j44j) about {k^,, M^) in a co-moving frame {X — c^^T, T) yields 

dT~k + kAc.{dx~k,dxM) = 

(45) 

dxM + d[<W^>]^{dxk,dxM)-c^dxM = 

14 



Performing a Fourier analysis, we look for solutions to equation ( I45p as 

{k,M) = e''^^+^'^\k,M) and exactly obtain the linear system (1330371) with 

A = Ao and (A;, M) = (/cq, Mq). As a result is also obtained the leading term 
in the expansion of the spectral curves associated to the stability of roll- 
waves, A(z/) = vk^K + 0(z/2). Thereby a necessary condition for roll-waves 
to be stable is that the Whitham's system (jUj) is hyperbolic. In order to be 
able to check this hyperbohcity in the limit 5 — )■ 0, we expand (l38i) . Recall 
from ^ that c = + 6M and < >= q + cM = + with c(°) 
and q^^^ independents of M. Thus the dispersion relation is 

A^ + {kj'^'ik^) - c(°)(A;.))A - (A;,c(°)(fc,)c(°)'(A;,) + 6MW'(A;J) = 

and system fH4|) is hyperbolic if 

A(A;,) := {c^'\k,) + kj^^\k,)f + 2Ak,qf\k,) > . 

We further carry out the spectral analysis of the second order Whitham's 
equations. The first of these two equations reads 

dTk + dx{kc{k,M)) = -edx[< k f >] (46) 

where f{k, M, drk, dxk, drM, dxM) is defined by < /, U' >= and 

A,m[/] = driU- <U>) + 3dx [U^- < >) 

+ 3 [k'^dlyyU + kdxk dp) + 5 {2kdxyU + dxk dyU) 
+ 5 {ied^xyyyU + ePdxk dp) . 

Notice that the e dependence disappears, e being rescaled to 1, by coming 
back to original variables {x,t) = {X/e,T/e). Obviously one could have 
equivalently considered u = 0{e) wavenumbers in the Fourier plane wave 
analysis of linearized equations. Afterwards, linearizing f l46p about (A;^, M^) 
yields in a co-moving frame equation 

dtk + k,dc,id,k,d,M) = -6k,d,<g> (47) 

where g = df^{k, M, dtk — c^^d^k, dxk, dtM — c^^d^M, d^M) satisfies both con- 
ditions < g, U[ >= and 

C4g] = dkUMk-cA~k) + dk{3U^-<3U^>)Jx~k 

+dM (3t/2- < 3f/2 >)JxM + 3 {klidkUj' + k^U'^ dxk 
+5 p: + 2K{duU.)' + Qkp':' + AklidkU,)'") dxk . 
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We search for plane waves solution in the form 0H'^)t+'^''*^ ^A;(z/), M(z/) j and 
expand A, k, M about z/ = as 

X{iy) = ukAo + {lyk^fXi + C»(z/3) , 
k{u) = ko + uk^ki + , 
M(z/) = Mo + vk^Mi + . 

Identifying 0{v) terms in fH7j) yields condition 

Ao/co + Mc^/co, Mo) = 

which is again equation fl33|l . We further identify O^u"^) terms and find 

AiA;o + AoA;i + A;,dc,(A;i, Mi) = (< /^^ > k^+ < > Mo) 

which is (l39l) . We focus now on the second Whitham's equation 

OtM + (9x < 3f/2 > = edx < QU{f- < / >) > -e6dxxM . (48) 

Rescaling fHHl) and linearizing about (/e^, M^) in a co-moving frame yield 

dtM + d[< 3t/' >].(5xfc, dxM) - c^dxM (49) 
= ax < 6U,ig- <g>)> -5dxxM . 

Identifying 0(z/) and C(z^^) terms yields respectively fl37j) and (H3l) . 

As 5 — )■ 0, the domain of validity of the Whitham's equations shrinks to 
and one has to consider modulation with wavenumbers whose magnitudes 
are of order 6. We analyze this situation in the next section. 

4 Stability of wave trains for S ^ 

In this section, we consider the stability of periodic wave trains as 5 — > 0. In 
what follows, we will write 6 = e6 where e is the characteristic modulation 
wavenumber. We first derive a set of three balance laws for modulated wave 
trains. In order to obtain a consistant spectral stability theory in the small 
wavenumber regime, and thus to be able to compare the Whitham's equations 
derived here and the spectral curves, we will have to consider wavenumbers 
u = 0{e). 
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4.1 Whitham's modulation equations 

Following section [3l we rescale (gKS) so that it reads 

dru + 6udxu + e^^^w + S {e'^dlu + e^dj^u) = . (50) 
We then search for modulated solutions in the form 

u(X,T) = uJ^e(l2ll,x,T]+euJ^e(^-.X.T]+0(.') , 



~fi{x,T) = ^„(x,r) + £43i(A',r) + 0(£') , 

with Uo{y;X,T) and Ui{y;X,T) 1-periodic in y. Inserting this ansatz into 
(150|) and identifying 0{e~^) yield after setting ko = dx^po and Qq = Ot^q 

^odyUQ + QkoUodyUo + k^dyUo = . (51) 

Thus Uq{- ] X,T) is a periodic traveling wave of (KdV). These solutions are 
given in term of elliptic functions (see section [2]), 

uo{y;X,T) = M^{X,T) + [Uo{k,{X,T),p,{X,T))]{y) , 
Co{X,T) = 6Mo{X,T) + do{ko{X,T),po{X,T)) , 

for some {Mo,po), with Uq and cq given by ffTOj) . From now on, we drop the 
subscript in the definitions flTU]) of U and c and keep it for the functions 
of paremeters {k,M,p) evaluated in {ko, Mq,pq). Moreover we will denot^ 
U{k,M,p) = M + jj{k,p) and c{k,M,p) = 6M + c{k,p). Compatibility 
condition drOxfo = dxdrfo yields then 

dTko + dxikociko,Mo,po)) = 0. (52) 

Next, we consider 0{e^) terms to obtain a system for {ko, Mo,po) in a closed 
form. Setting Qi = dr^pi and ki = dx4>i yields 

Co[ui] + QidyUo + ki [SkldyUo + GuodyUo) + OtUq + dx{^ul) 

+ 3ikldxyyUo + kodxko d^uo) + 6 {kld^uo + k^d^uo) = (53) 



with jC.k,p[f] = kdy i^{6ij{k,p) -d{k,p))fj + k^d^f. Equation ((53]) has a 

solution if and only two compatibility conditions are satisfied. Indeed C is 
Fredholm of index and the kernel of C^"^ is spanned by 1 and U (or U). 
The first compatibility condition yields 

Ot <uq> +dx < 3^0 > = 



^This notation is of course not consistant with the meaning of U and c up to now. 
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thus 

BtMo + 3 (9x < f^o > = . (54) 
The second compatibihty condition reads 

< StUq, Uo> + < dx{Sul),Uo > +3 < k^dxyyUo + k^dxh d^uo, Uq > 
+ 6{kl < d^uo, Uo > +k^, < 9>o, Uo >) = 0. 

where < ■ , • > is the usual scalar product on 1-periodic functions. Integrat- 
ing by parts yields equation 

<U',>+dx< 2Ul - hl{U'^^ > = 6kl < {Ul^r - kliUlir > . (55) 

As a result, one obtains one extra balance law in comparison to the case 
where 5 > is fixed. Indeed, in the case (5 — )■ 0, (gKS) is a perturbation of 
(KdV) and for 5 = 0, one exactly obtains the modulated system for (KdV), 
as it was obtained in [9]. However, in their derivation, they considered the 
energy equation 

dt (y ) + (2u' + udlu - = 

obtained from (KdV) by multiplying the equation by u, as an extra conserva- 
tion law and they averaged this conservation law to obtain fl55l) (when 5 = 0). 
Notice in particular that many conservation laws were available for such a 
process. Here, this extra equation seems to be obtained more naturally as a 
comptability condition in the process of computing asymptotic expansions of 
low frequency perturbations of the wave. The stationnary solutions to (l52l 
[54|55|l are exactly the parameters /cq^^o^Po such that the following relation 
holds 

< [u{ko,Po)')^ > = kl< {uiko,po)"y > . (56) 

It corresponds exactly to the solvability condition in order to compute an 
expansion of periodic solution to (gKS) as 5 — )■ 0. This restricts the set of 
periodic solutions to a two dimensional manifold just as in the case 6 > 
fixed. From the dynamical point of view, we recover a set of two equations 
by considering the relaxation limit 5 — )■ cxd. This yields the relaxed system 
fl52f54f56p which is a system of two conservation laws, just as in the case of 
fixed 5 > 0. 

Should conditions (15^ and (155|1 hold, one may solve equation 

jCo[vo]+dTUo + dx{Sul) 

+ 3{kldxyyUo + kodxko dyUo) + 6 {kldyUo + kgdyUo) = 
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with t>o, in a unique way when added condition^!!! 

<Vo,Ul^> = <Vo> = . 

Note that the fonction vq is odcF*l. Equation (jnSD reads then 

^o[ui -vq- hdkUo] + - hdk^o) Uq = . 

We denote Mi =< ui > and introduce pi the constant (with respect to y) 
such that Qi = dQQ[ki, Mi,pi]. Solutions to flS^ are then written as 

Ml = df/o[fci, Mi,pi] + t;o + 7if/o 

for some 71 independent of y. In what follows, we denote = dUo[ki, Mi,pi] 
and u° = Vq + 'JiUq, respectively even and odd parts of Ui. Compatability 
condition drOx't'i = dxdT(pi yields 

drh - dx {dno[kuMuPi]) = . (57) 

In order to obtain the full set of equations that governs the evolution of 
{ki,pi, Ml), one has to consider 0{e) terms in (150|) and finds 

Co[u2] + QidyU° + drul + dx{Quoul) 

+ Qdy{kiUoul + koulu") + 3{k^dxyyul + kodxko d^ul) 

+ 3 {2kokidxyyUo + dx{koki)dyyUo + klkid^u") (58) 

+ 6{kldyi + kld^ul + 2kokid^uo + 4:klkidluo) +F° = 

where F° is some odd function. First solvability condition for f l58p reads 

dr <ul> + dx < QuquI > = . 

Inserting ul into this latter equation yields 

dTMi + dx{d[<3U^>Uki,Mi,pi]) = 0. (59) 

The second solvability condition gives 

< OtuI, Uo> + < dxiQuoul), Uo> - < QkiUou° + Qk^ulul + Vtiul, Uq > 

- ^ (2 < kldxyul, U'^>+dx [kl) < dyul, U'^ >) 
-?>{dx< k^ki{U'^^ > -klki < dyul, >) (60) 

- 6 {kl < dyul, f/^ > -fco' < d'yul, Ui; >) 
-d{ki {2ko < (f/^)' > -Akl < {U^f >)) = . 

II See the description of the kernel of C in the next subsection. 
**0f course most of unobvious claims of oddness or evenness may be deduced from 
uniqueness. 
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Equations f lFr|59|60p govern the evolution of {ki, Mi,pi). 

If these equations are satisfied, we can further solve fISS]) and express U2 
in terms of (/c2,M2,p2) and {ki,pi, Mi,Ui)i=o^i. Actually, by induction, one 
can proceed with expansions up to any order with respect to e. 

Unlike the case 5 > fixed, by summing fl52ll54f55l) and £ fl57ll59f60p 
we can not obtain a closed system governing the evolution of {k,M,p) = 
{ko, Mq,pq) + e{ki, Mi,pi) up to order O^e"^) terms. The first two equations 
would indeed close into 

dxk - dx^ = , 
OtM + (9x < 3f/2 > = , 

but the last equation would not. However, we will relate the stability of 
periodic wave trains in the limit 5 — > with a wave analysis of the two first 
steps of the modulation. 

4.2 Bloch analysis of spectral stability 

In this section, we analyse the spectral stability of periodic wave trains in 
the limit 5 — )■ 0. Recall that they expand as 

U\y;k,M) = M + U{y;k,p) + 6U,{y;k,p) + 6^U2{y;k,p) + 0{6') 
c\k,M) = 6M + c{k,p) + 6^C2{k,p) + 0{6^) 

with parameter p obtained from k through equation f lT2|) . Moreover stand 
the following relations 

eif!' + k({6u- ~c)Ui)' + eu" + k'^u"" = , 

) . J . . . (61) 

k^V'^ + k f (6f/ - c)[/2 + Wl - C2f/j + k^V'i + A;Vf' = . 

Notice U\ is odd and [/2 is even. 

The spectral problem associated to f lSOp is written in a co-moving frame 

Kv)u^K [I + v) ((6f/f - cf) u) + kl {± + vfu ^^^^ 
+ ^''^{kl{± + ufu + kt{± + uYu) = 

with a 1-periodic u. In the limit e — )■ the analysis of f l62p has been carried 
out in [1] both for v = 0{1) and for u = 0{e) (in an appendix). Here we 
recall only the latter one since this is the one to be interpreted in terms 
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of linearized modulation equations. Thus we denote z/ = z/e and expand 
eigenfunctions and eigenvalues as 

A(z/) = eVXo + {eVyXi + 0{e^) , u = Uo + eui + e^U2 + 0{e'^) . 

Identifying 0{£^) terms in flB^ . one finds 

L[u^] = A;, ((6t/.-cjMo)' + A;X" = 0. 
The two dimensional kernel of is spanned by U[ and 

OpC^ 

Thus Mo = kiU^ + MoKt for some (A;i,Mo). Notice < Uq >= Mq. Next, we 
consider 0{e) terms 

-C^iii] = kj (gUu Uo)' +s iklii'^ + ktii'^') 

_ q — (~ ~ ^ \ ^^"^^ 

+ "ivklxi!^ + vk^ ( AqMo + (6t4 — Cj,)'Uo ) 

where again Aq stands for Xo/ki,. Solvability conditions are < fl63|) r h » >= 
and < [/^. fl63|) r h » >= 0. They reduce to Mq = 0. Equation fl63l) then reads 



- uk^kidkUi, - ki5U[ ,\ + z/fc^fci ^Aq + KdkC^ U'^ = Q . 
We introduce Mi =< -ui > then pi the constant such that 

Ao/ci + Mc.(A;i,Mi,pi) = (64) 
where c(A;, M, = 6M + c{k,p). Therefore accordingly 

Ml = Vk,dU,{k,,Mi,pi) + k,6Ul, + k2Ul (65) 

where again U{k,M,p) = M + U{k,p). To determine (Aq, fci, Mi,pi), one 
must look at fl62|) up to order 0{e) and then finds 



-A[m2] = i^iXoUi + uXiUo) + kJ (eUi^^ui + (5(6f/2,^ - C2,^)mo^ 

+ 3A;3l7 (m'Z + Vu'q) + z/A;, ((6f/, - c.)mi + 65 t/i,.Mo) ^^^^ 
+ Ski {u'l + 2z7m'o) + Sktiu'l" + AVu'^) . 
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Solvability conditions are < 



r.h.s. 



>= and < 



into account fl^Sjl first condition reads 

(Ao-c,)Mi + d[<3f/2>],(fci,Mi,pi) = 0. 
With odd and even parts notations, equation ( l66|) is also written 



h.s. >= 0. Taking 
(67) 



where L° is some odd term. Using the second equation of fl6Tl) and the fact 
that 

< f/., k^il!^", + A;. f(6f/. - c.)?7L)' > = 



(coming from £*f/^ = 0), the last solvability condition is thus written 

= 77A;, [(Ao - cjd < if/2 >. +d < 2f/3 - |P(f/')2 >,] {k,,M^,p,) 
-5 d < A;2(f/')2 - A;4([/")2 >^ (k^,M^,p^) 
-QkJ <U,,,UidU,ik,,Mi,pi)> 

+ M f (Ao - < t7{ ,f/. > - < 3klU'^^ + 6Ui,M., U: > 



(68) 



As a result, fl64|67|68p forms a full spectral problem for eigenvalue Aq 
and eigenvector {ki, Mi,pi). Therefore Aq must satisfy a dispersion relation 
E{Xq) = with E a polynom of degree 3. This dispersion is analogous to the 
one found in [Ij. In the following, we interpret the system fl64f67f68p with 
the help of a wave analysis of modulation systems derived in the previous 
section. 

4.3 Spectral validity of Whitham's equations 

Let us first remark that a byproduct of the previous analysis is the direct 
spectral validation of modulated Whitham's equations for (KdV). We recover 
(KdV) setting 5 = in previous computations. The linear spectral problem 
fl64f67f68p reads then 

I (Ao - c^)dA;^ + d[A;c],, (A;o,Mo,po) = 

(Ao-c,)dM, + d[<3f/2 >]J(A;o,Mo,po) = 



^^(Ao-c.)d[< f/2 >]. + d[<2U' 



> 



\k < 



([/')'>] J(A;o, Mo, po) 
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whereas the modulated system for (KdV), obtained by setting 6 
is 



drko + dxikoCo) = 



in 



dxMo + dx < Sf/n^ >= 



(69) 



dT<U^> + dx< 2f/3 - |A;2(?7^)2 >= o 



L 2 



There is no difficulty here to see that the spectral problem described above 
and the one obtain when considering the stability of steady solutions to ( 1^ 
are the same. This validates modulated equations for (KdV) on a spectral 
level. Moreover, hyperbolicity of modulated equations fj69|) is thus a neces- 
sary condition for periodic travelling waves solutions to (KdV) to be stable. 

The situation is trickier for (gKS): from the spectral analysis carried 
out in the previous section, we see that the correction Ui to U plays a sig- 
nificant role. As a result, the linearization of f lS^54|55l) at a steady solu- 
tion is clearly not sufficient and we need to take into account higher order 
modulations. Actually the right way to search for "solutions" in the form 
(fc*, Mj,,p^) +e(/ci, Mi,pi), with {k-^, M^,p^) constant and {ki, Mi,pi) a plane 
wave, is to set in the modulation process (/cq, Mq,pq) = {k^, M^,p^) and look 
for {ki, Mi,pi) as a plane wave. Obviously, when 5 = 0, it coindides with 
the direct approach. 

Setting (fco, Mo,pq) = (A;^, M^,p^,) (with p^ obtained through f|T2|) ). equa- 
tions (1571) and fl59l) turn into 



drk^ + d[A;c]A(fci,Mi,pi) = 
OtMi + d<3U^ >^dx{ki,Mi,pi) = 



(70) 
(71) 



and, since then vq = using 

<[/„£,(df/:(A;i,Mi,pi))>= 0, 

condition ( 160|1 reads 

= d<^U^>, ar(fci,Mi,pi) + d < 2f/3 - |A;2(f/')2 9x(fci,Mi,pi) 
-5d< k\U'f - k\U"f >. 9x(A;i,Mi,pi) 
-6M < t/i,.f/:d[A;f/].(A;i,Mi,pi) > 
- M < 3A;2f/(;^, U'_>-'6< Ul^U^ > d[kc]^{ki, Mi,pi) . 

Searching, in a co-moving frame, for plane- wave solutions to 
leads exactly to the spectral problem f l6^67f68p . This validates at the spec- 
tral level the modulation equations in the limit (5 — )■ 0. 



(72) 
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5 Conclusion 



In this paper, we derived Whitham's equations for the modulations of pe- 
riodic travehng waves, solutions to generalized Kuramoto Sivashinski (gKS) 
equations. We proved that this set of equations gives necessary conditions 
of first and second order for the spectral stability of travelling waves in the 
low frequency regime. Indeed, the hyperbolicity of Whitham's first order 
equations implies the tangency of spectral curves at the origin whereas the 
parabolicity of the second order Whitham's system ensures that the spectral 
curves lie in the stable spectral region. We also reinterpreted the spectral 
analysis carried out in pQ in the (KdV) limit (5—7-0 with the help of a set 
of Whitham's equations. As a byproduct, we have obtained a set Whitam's 
equations for modulated wavetrains in (KdV) equation in a more natural 
way than in where an additional equation for the energy was introduced 
to close the Whitham's set of equations. 

It is worth noting that, following the approach already introduced in [B] in 
the case of reaction diffusion equations, one could have also derived a viscous 
Burgers equation by choosing the scaling: 

u{X, T) = e'u,{^^^, X - c,T, eT),), 
0(X,T) = MX,T) + J2e''<p\X - CgT,eT), 

k>l 

where 0^, is the phase associated to a fixed periodic wave train Ui,{.,k^,, M^) 
and Cg is one of the characteristic speeds associated to the inviscid Whitham's 
system. This is the first step towards the construction of approximate solu- 
tions to (gKS) on asymptotically large time interval 0{T*/e'^) with arbitrary 
T* > 0. Furthermore, as viscous Burgers equations have shock structures, 
we expect that we will be able to construct solutions of shock type in (gKS) 
in the form of a travelling wave relating two families of periodic wave train 
at — oo and +oo. This approach will be developped in a forthcoming paper. 

Another problem of interest is the stability of periodic travelling waves, solu- 
tions to (gKS): we expect that, using the approach introduced in [8], one can 
prove the nonlinear stability of periodic wavetrains solutions to (gKS) under 
suitable spectral assumptions. Concerning this spectral stability issue, there 
should be a range of periods for which periodic wavetrains are stable, just as 
in the case of roll-waves solutions to shallow water equations. Indeed, at the 
onset of periodic wave train (close to the Hopf bifurcation), steady solutions 
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are unstable and by perturbation arguments, we shall prove the same result 
for small amplitude solutions. On the other hand, solitary waves are also 
unstable since, for (gKS), constants are unstable under low frequency per- 
turbations. This heuristic approach was partially validated in [1] in the limit 
5 — )■ but the authors only considered perturbations with zero mean. This 
analysis should be completed for general perturbations in order to conclude 
to the full stability of periodic wavetrains. 
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